Phase slips and phase synchronization of coupled oscillators 
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The behaviors of coupled oscillators, each of which has periodic motion with random natural 
frequency in the absence of coupling, are investigated. Some novel collective phenomena are revealed. 
At the onset of instability of the phase-locking state, simultaneous phase slips of all oscillators and 
quantized phase shifts in these phase slips are observed. By increasing the coupling, a bifurcation 
tree from high-dimensional quasiperiodicity to chaos to quasiperiodicity and periodicity is found. 
Different orders of phase synchronizations of chaotic oscillators and chaotic clusters play the key 
role for constructing this tree structure. 
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The investigation of coupled oscillators has attracted 
constant interest for many decades [1] . The rich collective 
behaviors of these systems, such as mutual entrainment, 
self-synchronization, and so on, are observed in many 
fields, e.g., coupled laser systems, Josephson junction ar- 
rays, biological and chemical oscillators etc. [2-10]. In 
early studies, interest was focused on coupled oscillators, 
of which each is periodic without coupling. Recently, the 
investigation has been extended to coupled chaotic sys- 
tems (i.e., individual systems are chaotic without cou- 
pling). Significant phenomena were found, such as phase 
synchronization of two mutually coupled chaotic oscilla- 
tors [11, 12], and clustering and cluster-cluster synchro- 
nization of multiple coupled chaotic units for local [13] 
and global [14] couplings. 

In this letter, we study the following N coupled oscil- 
lators with the nearest coupling, 
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i = 1,2,. ..,N, where K, 9i and Ui are the coupling 
strength, the angle of modulo 27r and the natural fre- 
quency of the z-th oscillator, respectively. Model (1) has 
been extensively investigated in the past several decades. 
Here we concentrate on the dynamical behavior of the 
system. In particular, we are interested in the charac- 
teristic features of the motions of individual oscillators, 
i.e., the microscopic motions, in the regime of desyn- 
chronization of the phase-locking state, which have not 
yet been well investigated by the previous works. Sev- 
eral novel features of this system are found. First, we 
find simultaneous phase slips of all oscillators at the 
onset of desynchronization of the phase locking state, 
and quantized phase shifts in these slips are observed. 
Moreover, we find the interesting cascade behavior of 
coupling-induced chaos and a nice tree structure of tran- 
sitions from qausiperiodicity to chaos to qausiperiodicity 
and periodicity. Then rich behaviors of synchronizations 
between chaotic oscillators and chaotic clusters, which 
have attracted much attention recently for the coupled 



Rossler and Lorenz oscillators, can be also identified for 
the rather old as well as popular system (1), of which the 
individual oscillators have simple periodic motions with- 
out coupling. These findings greatly enlarge the appli- 
cation perspectives of chaos synchronization. These fea- 
tures of phase dynamics are expected to be observable in 
practical systems by experiments, such as coupled laser 
arrays, Josephson junction chains and coupled electrical 
circuits. 

In Eq. (1) the periodic boundary condition 0i + N(t) = 
0i(t) is applied. Without losing generality we scale tOi 
such that 
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It is well known that for a given N and {oji, i = 1, ...N}, 
there is a critical coupling K = K c . For K > K Cl phase- 
locking can be observed, and then we have {9i= 0, i = 
1,...,N} , and each di is locked to a fixed value. For 
K < K c , no phase-locking exists, and 0i it) are no longer 
zero. In [8-10], it is found that if we define an average 
frequency as 
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synchronization between different oscillators, in the sense 
of u>i=u)j, i 7^ j, can be observed in the region where 
strict phase- locking of 0i= is broken. It is interest- 
ing to investigate how the various oscillators are led to 
complete synchronization (phase-locking for K > K L ) via 
a sequence of bifurcations by increasing the coupling K 
from K = 0. 

In order to get a general idea about the global behavior 
of the system, we first measure the following two positive 
macroscopic quantities R and fi: 
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It is clear that R is time-dependent beyond the phase- 
locking region. Then we further measure its time average 
as an order parameter: 

(R)= Hm i [ T R(t)dt, (5) 

In Figs. 1(a) and (c), we plot Q vs. K for N — 5 and 15, 
respectively. In both cases, natural frequencies are ran- 
domly chosen from a normal Gaussian distribution. The 
actual frequencies can be seen in Fig. 3(a) for N = 5 
and Fig. 3(b) for N — 15 at K = 0, respectively. These 
natural frequencies are used for all figures in this letter. 
We found, for the given natural frequencies, K c = 5.08 
for N = 5 and K c = 6.22 for N = 15. When K > K c 
we have identically = 0, and then complete synchro- 
nization (phase- locking) is justified. In Figs. 1(b) and (d) 
we do the same as (a) and (c), respectively, with the 
measured quantity replaced by < R >. In computing 
Figs.l, initial conditions of &j(0) are randomly chosen, 
and then in the shaded region of Fig. 1(d) the coex- 
istence of multiple attractors of phase-locking states is 
justified for K > K c . In Fig.l, the quantity T in Eqs. 
(4) and (5) is taken sufficiently long in our simulations 
so that the fluctuations due to finite T are invisible. It 
is interesting to observe that several discontinuities ap- 
pear in the region K < K Cl indicating that, apart from 
the apparent phase-locking transition, some additional 
transitions exist even before K c . 

For the microscopic quantities, it is natural to study 
the velocities of various oscillators. In Figs. 2(a) and 
(b), we present the motions of Oi vs. t for N = 5 and 
different AT's. For K = 0, 6i must be equal to the con- 
stant tUi, and for small K [see (a)] 0i varies oscillatorily 
around its natural frequency. As K increases, the oscilla- 
tions of Oi become large, and the oscillation centers of all 
oscillators shift close to each other. An interesting fea- 
ture is that near the onset of synchronization [see (b)], 
we can find simultaneous phase slips, i.e., all oscillators 
keep in a the phase-locking condition (OFF-state) for a 
long time, and then simultaneous bursts of all oscillators 
(ON-state) break the locking state. After a short firings 
all oscillators calm down again simultaneously to another 
phase-locking state, and then repeat the same process pe- 
riodically. As K gets closer to K c , the length of the OFF 
state r becomes longer. We find a clear scaling between 
r and K c — K : 

toz(K c -K)-°- 5 . (6) 

The above features of synchronized actions of phase 
slips can be well understood by an intuitive explana- 
tion. Suppose various oscillators can be locked to a set of 
phases Oi (K) for K > K c . Then all solutions satisfying 

Bi+i {K,m}- ~Qi (K,m) =Q l+ i (K)- Oi (K) + 2irrrn 

(7) 



with m = (mi, ...,m,j, ...mjv) and rrii being any integer, 
must also be phase- locking solutions of (1). By reducing 
K lower than K c , all the above solutions lose their sta- 
bility via the saddle-node bifurcation. At K = K c , there 
exists a heteroclinic path linking some of the above so- 
lutions, which has the lowest potential, and is attracting 
(the existence of such a heteroclinic path can be rigor- 
ously proven for N = 2). For K < K c , and \K — K c \ < 
1, the system takes such a periodic path, which stays in 
the vicinity of one of the above stationary solution for a 
long time, and escapes away from this solution (simulta- 
neous phase slips for all oscillators), and then approaches 
to the vicinity of the next stationary solution alone the 
heteroclinic path of AT = K c : this produces the peri- 
odic pulses of Fig. 2(b). The scaling property of time 
length of the OFF-state r can be also computed since 
for the saddle-node bifurcation we have a universal form 
x= (K c — K) + x 2 . The time r for x to move from x = 
to x -» oo reads r oc J °° {Kc _% +x 2 = 2 ^^_ K \ this 
explains our observation of the scaling law of Eq. (6). 

It is interesting to compute the phase shifts for various 
oscillators in each phase slip pulse in Figs. 2(b). Let A6i 
represent the phase shift of Oi during each pulse. From 
(1) and (2) we have 
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We argue that any two adjacent fixed points in a hete- 
roclinic path at K = K c take m* = or ±1 in Eqs. (7), 
i.e., 

A0 l+ i - AOi = or ± 2tt (9) 

for i = 1, N — 1. Considering both conditions (8) and 
(9), AOi can take only quantized values 

A, = 0,4.4,..., ± ^, ±2 , (10, 

The concrete value for each AOi depends on the partic- 
ular distribution of Wj. In Figs. 2(c) and (d) the phase 
shifts observed in all pulses fully agree with our heuris- 
tic argument. Further considering Fig. 3, the phase shifts 
in Fig. 2 can be exactly predicted (see [15]). One more 
conclusion from the above analysis is that at the onset 
of instability of the phase-locking state, all velocities u)j 
have the scaling Wjoc \/ K c — K, and the ratios between 
different u)j must be quantized to discrete rational num- 
bers. 

As we further reduce K to values considerably smaller 
than K c , no more phase- locking exists, and no appar- 
ent synchronization can be observed directly for Oj (t). 
However, some other implicit synchronization — phase 
synchronization — can be still found. To get a general 
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idea, we plot <Dj defined in (3) vs. K for N — 5 and 
15 in Fig. 3(a) and (b), respectively, by varying K from 
K = to K > K c . Interesting behavior of transition 
tree for phase synchronization is clearly shown. Though 
some synchronizations can be expected from the obser- 
vations of frequency plateaus reported in previous pa- 
pers [8-10], a number of characteristic features revealed 
in these trees are novel and interesting. Three kinds 
of transitions can be observed in these trees. First, if 
two adjacent oscillators (or adjacent clusters of oscilla- 
tors) have close frequencies, they can be easily synchro- 
nized by increasing K. In this case, one always finds two 
branches merging to a single one (indicated by A). Sec- 
ond, if two non-adjacent oscillators (or two non-adjacent 
clusters) have close frequencies while the oscillators be- 
tween them have considerably different frequencies, one 
can find the non- adjacent oscillators can be also synchro- 
nized to each other, i.e., non-local clusters can be formed, 
and these non-local clusters can quickly bring the oscilla- 
tors between them to the synchronized status, and form a 
solidly larger synchronized cluster. In this case, the tran- 
sition may be from two to one or from multiple branches 
to one (indicated by B). An oscillator, which is synchro- 
nized to a cluster for certain K, can be desynchronized 
from the original cluster by increasing K. This desyn- 
chronization always happens at an edge oscillator of a 
cluster, due to the competition between two neighbor 
clusters (indicated by C, see 2nd and 3rd oscillators of 
Fig. 3(b)). It is obvious that C is the inverse of A. The 
transitions of types B and C have never been realized 
before. 

The most interesting and novel finding is the nature 
of these synchronizations. The motions in these synchro- 
nization trees may be very different. They can be peri- 
odic, quasiperiodic and chaotic. In Fig. 4(a), we take 
N = 15 and plot the largest Lyapunov exponent of the 
system vs. K . In a large interval of K, we find the posi- 
tive Lyapunov exponent, indicating chaos. Thus, in this 
region phase synchronizations of chaotic oscillators are 
identified. Recently, the phase synchronization of cou- 
pled chaotic systems has attracted great attention [11]. 
Here the major difference between our system and the 
previous chaos synchronization is that our oscillators are 
periodic without coupling, and chaos is induced by the 
coupling, while in the latter case the individual systems 
are chaotic without coupling. Moreover, we find that 
these coupling-induced chaotic motions of different os- 
cillators may have different levels of synchronizations by 
varying the coupling, and the cascade of synchronization 
forms a tree-like structure of Fig. 3. 

In Figs. 4(b) to (f) we plot the maps of 6i (n) to 
01 (n + 1), where Oi {n) is the Oi value at the time t 
when 0\(t) crosses the angles 2nw with n being an inte- 
ger. For K > K Cl we have fixed point solution, and the 
map is fixed at Q\ (n) =0i (n + 1) = 0. For K slightly 



smaller than K c we have periodic solution represented 
by the finite number of dots in Fig. 4(b). The period 
8 can be easily understood from Fig. 2(d). The period 
of the total system is 15t, where r is the time interval 
between the two adjacent slips, and the change of of 9\ 
in 15r is — 167r. This leads to the period-8 solution of 
Fig. 4(b). Two- frequency torus can be identified in the 
three-cluster regime [see Fig. 4(c)]. For very small K, 
we can find high-dimensional quasiperiodicity [e.g., Fig. 
4(f)]. Chaos is prevailing in the K region between Fig. 
4(c) and 4(f) [see Fig. 4(d) and (e), and Fig. 4(a)]. 
Many periodic windows are found in the quasiperiodic 
and chaotic regions, which will be investigated in detail 
in our forthcoming extended paper. 

The entire variation from high-dimensional quasiperi- 
odicity (for very small K) to periodic motion (K < K c , 
\K — K c \ <C 1) through various orders of chaos synchro- 
nizations can be vitally seen in Figs. 3 and 4. Starting 
from the high-dimensional quasiperiodicity for K <C 1, 
by increasing K various neighboring oscillators with close 
frequencies start to form clusters via phase synchroniza- 
tion, and chaos is induced near the first synchronization. 
Then in each cluster, different oscillators perform dif- 
ferent chaotic motions, while having identical winding 
number. The winding numbers for different chaotic clus- 
ters are different. When further increasing K , adjacent 
chaotic clusters can be synchronized to form larger clus- 
ters, until two large clusters are formed when the motion 
becomes periodic. This tree picture of the transitions is 
expected to be the same for general large number of cou- 
pled non- identical oscillators, which are periodic in bare 
case. 

Since system (1) qualitatively describes practical sit- 
uations in wide fields, ranging from physics, chemistry 
to biology, the findings in this letter are expected to be 
of general significance, and they can be used for under- 
standing the mechanisms of rich collective behaviors of 
coupled systems. Laser arrays, Josephson junction chains 
and coupled electrical circuits may be ideal candidates 
for experimentally revealing the features explored in this 
letter. 
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From Fig.3(a) and Eqs.(7) and (8), we have A6> * = 

1,2,4,5 

-A<9 3 = 4A&, i = l,2, 4, 5, and A6 2 - A8 3 = 2tt, which 
lead to A6 3 = and A6>; = ^, i = 1,2,4,5. The 

same computation from Fig. 3(b) leads to 8A#, 
7A6W2-3, A6> 4 - A6> 3 = 2tt, producing A6>; = 



i = 4 - 11 and A6„ = 
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FIG. 1. Q. and < R > defined in Eqs. (4) and (5) plotted 
vs. K for N = 5 and 15. The natural frequencies are shown 
in Fig. 3 at K = 0. K c = 5.08 and 6.22 for N = 15 and 
N = 15. For K > K c , phase-locking can be observed, and 
the shaded region in (d) indicates the coexistence of multi- 
ple attractors of phase-locking states. Initial conditions are 
randomly chosen for each K detected. 
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FIG. 2. (a), (b): The time evolutions of Oi plotted for 
N — 5 and two different K. Near the onset of the instability 
of the phase- locking state [(b)], simultaneous phase slips of 
all oscillators are clearly seen, (c), (d): Some 8i(t) plotted 
against t near K c . (c) K = 5.07, N = 5; (d) K = 6.212, 
iV = 15. Quantization of phase shifts for all oscillators at 
each phase slip and the pulse are justified. 
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FIG. 3. Transition trees of synchronization for averaged 
velocities versus K. (a) N = 5, (b) N = 15. Note the exis- 
tence of three kinds of transitions labeled by A, B and C. 
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FIG. 4. (a) The largest Lyapunov exponent of the cou- 
pled oscillators plotted against K, N — 15. (b)-(f) give the 
maps of 9i (n + 1) vs. Bi (n) where 6i (n) is the value of 
0i (t) when 6<i crosses 2im. (b) K = 6.0, the motion is pe- 
riodic, (c) K — 3.0, the motion is quasiperiodic with two 
irreducible frequencies, (d), (e): K = 1.5 and 0.9, respec- 
tively. The motions are chaotic, (f) K — 0.01, the motion is 
high-dimensional quasiperiodicity. 
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